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I.  SUMMARY  OF  TECHNICAL  ACCOMPLISHMENTS 


The  purpose  of  the  research  under  this  contract  was  to  study  the  processes  by  which  sto¬ 
chastic  motion  is  self-generated  in  deterministic  systems,  and  the  consequences  of  the  resulting 
stochasticity.  We  consider  representative  problems,  each  of  which  incorporates  one  or  more  of 
the  main  underlying  physical  phenomena.  Two  classes  of  phenomena  that  we  have  studied  under 
this  contract  are  diffusion  across  overlapping  resonances,  and  diffusion  along  resonances,  in  more 
than  two  degrees  of  freedom  (Arnold  and  modulational  diffusion).  The  former  has  been  studied  in 
considerable  depth,  and  new  insights  involve  the  fine  structure  of  the  phase  space  or  more  exact 
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calculations.  The  latter  phenomena  still  have  many  poorly  understood  aspects. 

1.  Hamiltonian  Renormalization  Theory  • 

i  * 

A  systematic  Hamiltonian  renormalization  method  using  regular  (secular)  perturbation 
theory  has  been  formulated  to  determine  at  what  energy  the  resonance  structure  of  the  system 
repeats  itself  to  higher  orders  [1].  A  self-similar  property  of  Hamiltonian  maps  has  been  deter¬ 
mined  for  p  iterations  of  the  mapping,  where  p  is  an  arbitrary  integer.  The  procedure  involves  the 
order-by-order  renormalization  of  resonances,  in  which  each  resonance,  averaged  over  the  fast 
variables,  locally  approximates  the  motion  of  a  pendulum.  A  standard  mapping  is  reconstituted 
for  each  resonance.  Locally,  this  becomes  the  approximate  description  of  any  mapping  near  a 
resonance.  The  value  of  the  stochasticity  parameter  K ,  for  self-similar  behavior  about  a  p- 
interaction  fixed  point,  is  then  determined. 

2.  Corrections  to  Quasllinear  Theory 
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Higher  order  corrections  to  the  quasilinear  diffusion  coefficient  can  be  obtained  for  the  stan¬ 
dard  (Chirikov -Taylor)  map.  We  have  studied  whether  these  corrections  can  be  applied  to  other 
maps  which  can  be  locally  approximated  by  the  standard  map.  We  have  analytically  obtained  a 
local  diffusion  coefficient  as  a  function  of  action  for  the  Fermi  map  [21  which  displays  the  charac¬ 
teristic  oscillations  in  action.  We  have  used  this  diffusion  coefficient  in  the  Fokker-Planck  equa- 
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tion  to  obtain  numerically  the  diffusion  of  an  initial  delta  function  distribution  in  action.  We  have 
compared  this  with  a  numerical  calculation  of  the  diffusion  for  the  same  initial  distribution 
obtained  by  direct  iteration  of  the  mapping  equations.  The  comparison  demonstrates  that  the  local 
approximation  can  be  used  to  obtain  the  higher  order  corrections  to  the  diffusion  over  much  of 
the  stochastic  region  for  the  Fermi  map  [2J. 

For  some  values  of  parameters  the  standard  mapping  exhibits  accelerator  modes  for  which 
the  diffusion  becomes  mfinite.  For  the  more  generic  mappings  the  modes  become  quasi-modes 
and  the  diffusion  is  enhanced  over  values  of  the  action  for  which  the  quasi-modes  are  present.  We 
have  studied  this  phenomenon  and  developed  a  procedure  for  calculating  the  diffusion  in  the 
presence  of  the  quasi-accelerator  modes  by  treating  them  as  sources  and  sinks  in  the  Fokker- 
Planck  equation  [3]. 

We  have  completed  a  study  of  stochasticity  and  resonances  in  the  two  beam  accelerator  [4], 
The  two  beam  accelerator  is  a  proposed  linear  e+e~  machine  employing  2  x  500  tapered  free  elec¬ 
tron  lasers  (FELs)  as  a  source  of  microwave  power.  A  high  current,  low  energy  (20  MeV)  beam 
drives  the  FELs,  producing  1  cm  microwaves.  These  microwaves  are  used  to  accelerate  a  low 
current  beam  to  energies  of  1  TeV.  High  efficiencies  are  obtained  by  re-accelerating  the  low 
energy  beam  after  passing  through  each  FEL.  Restoring  the  2  MeV/particle  lost  on  passing 
through  each  FEL  section  (using  an  induction  accelerator)  avoids  wasting  the  18  MeV/particle 
left  in  the  beam.  The  periodic  nature  of  the  acceleration  and  deceleration  of  the  low  energy  beam 
can  lead  to  stochastic  motion.  We  have  developed  design  criteria  that  ensure  that  detrapping  of 
the  low  energy  beam  by  this  stochasticity  is  minimal,  preserving  the  high  efficiencies  inherent  in 
a  single  tapered  FEL.  Numerical  integrations  of  a  one  dimensional  model  for  the  FELs  were  per¬ 
formed,  showing  various  degrees  of  detrapping.  We  also  explored  the  effects  of  islands  produced 
by  resonances  between  the  periodic  acceleration  and  the  trapped  particle  motion.  These  islands 
represent  coherent  motion  of  large  numbers  of  trapped  particles,  leading  to  oscillations  in  the 
power  output  of  the  FELs. 
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3.  Arnold  and  Modulational  Diffusion 

Our  studies  of  Arnold  diffusion  in  the  billiards  problem  have  been  completed  and  a  review 
has  appeared  [5].  A  major  advance  has  been  the  application  to  the  ECRH  problem.  For  one 
specific  case  of  three  strong  resonances  which  are  approximately  equally  spaced  in  frequency,  the 
thin  layer  Arnold  diffusion  along  a  resonance  layer  has  been  treated  theoretically,  yielding  diffu¬ 
sion  coefficients  in  reasonable  agreement  with  numerically  determined  values  [6],  This  mechan¬ 
ism,  which  can  cause  particles  to  diffusively  heat  across  an  adiabatic  barrier,  is  found  to  be  some¬ 
what  weaker  than  collisional  diffusion  for  typical  plasma  confinement  devices.  However,  Arnold 
diffusion  could  prove  significant  for  additional  heating  of  plasmas  having  higher  temperature  (or 
lower  density)  [7], 

In  modulational  diffusion  an  oscillation  in  one  degree  of  freedom  is  modulated  at  a  slow 
driving  frequency,  yielding  a  broad  band  of  stochasticity.  The  coupling  of  this  stochastic  motion 
to  other  degrees  of  freedom  is  known  as  modulational  diffusion.  The  diffusion  coefficient  has 
been  evaluated  analytically,  without  the  need  for  adjustable  fitting  parameters.  The  coefficient  has 
been  compared  with  numerical  calculations  over  a  wide  range  of  system  parameters  and  for  diffu¬ 
sion  rates  varying  over  many  orders  of  magnitude.  The  numerical  and  analytical  results  are  in 
good  agreement  [8], 

4.  Intrinsically  Stochastic  Heating  Mechanisms 

A  study  of  two  frequency  heating  in  the  Fermi  map  showed  that  using  two  frequencies  gave 
an  approximately  twofold  increase  in  energy  for  the  position  of  the  lowest  KAM  barrier  to  sto¬ 
chastic  heating  when  the  resonances  associated  with  the  two  frequencies  were  properly  interlaced. 
We  have  extended  this  work  to  consider  the  more  practical  problem  of  two-frequency  electron 
cyclotron  resonance  heating  (ECRH)  in  mirror-confined  plasmas.  The  ECRH  is  modeled  by 
means  of  a  four  dimensional  simplectic  mapping  denved  from  the  nonrelativistic  single  particle 
energy  change  due  to  the  spatially  separate  resonance  zones.  Fixed  points  have  been  located  and 
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their  linear  stability  determined  analytically.  Resonances  in  action  space  have  been  calculated 
and  used  to  obtain  the  adiabatic  barrier.  Analytical  expressions  were  derived  for  the  quasilinear 
diffusion  coefficients  in  the  stochastic  regime  below  the  adiabatic  barrier  and  found  to  agree  well 
with  numerical  calculations.  Intrinsic  diffusion  in  the  parallel  energy  is  shown  to  lead  to  axial 
losses  that  can  rival  those  induced  by  collisions  alone.  Detailed  papers  describing  these  results 
have  been  published  [6,7],  A  new  phase  space  reconnection  mechanism  was  uncovered  and 
explored  in  a  simpler  system  [9], 
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